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PREFACE 



The object of this study is to determine the coefficient of eddy 
conduction, its variation with height, and its diurnal variation in the 
surface layer. 

This paper -.ms prepared at the U, S, llaval Postgraduate School, 
Monterey, California during the period December 1?4? to May 1^50, 

The author expresses his appreciation and thanks to Associate 
Professor Frank L. Martin, who gave valuable advice and guidance in 
the development of this paper. 
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I. INTRODUCTION 



During I748-I749 sevoral articles have been published by the 
University of Texas giving half-hourly observations of both tempera- 
ture and wind at various levels from the ground up to and including 
288 feet for periods of approximately one day* This suggested tho 
problera of determining tho coefficient of eddy conductivity using the 
observed lapse rate data; and, the coefficient of viscosity by using 
the observed vrind profile data. These tvro coefficients were found, in 
by Sverdrup Cn] to be equal in tho atmosphere, but other eddy 
diffusion coefficients wero found by him to be unequal. On the other 
hand, Pctterssen and Swinbank [10] , in 1947, suggest that the former 
coefficient is more nearly related to the latter by the factor .65 in 
the free atmosphere. The question appears to require another critical 
examination; and, with the availability of new and fairly accurate data, 
it seemed to be a problem well worthy of study. In actuality, the problem 
of determining the diurnal variation of the coefficient of eddy conductivity, 
K, near the ground proved to be so complex that this was all that time per- 
mitted in this dissertation. It will be necessary to answer the question 
of the equality of the two coefficients at some later date. However, tho 
determination of values of the coefficient of eddy conductivity and its 
diurnal variation, determined purely from heat transfer considerations, 
is of some practical interest as a check on previously determined values. 

For example, Stewart |]lj indicates a possible range of 1CP cm 2 sec"l to 
102 cm^ sec~3- for K within the surface layer, depending upon the stability 
of the layer. 
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This study gives explicit values of K, with respect to local civil 
time, corresponding to a particular synoptic situation, that for the 
period 27-28 September I94O in the vicinity of l.'anor , Texas, Another 
important result of this study proved to be the determination of ratos 
of radiative heating and cooling of the layers of air near the ground at 
various times covering a full day. This yielded some surprisingly large 
values. These results could be used to solve a problem proposed by 
Brunt DC concerning the nagnitudo of the coefficient, I",,, of radiative 
diffusivity. Brunt originally computed to be about 650 cm^ soc -l but 
in 1?40, using the new absorption coefficients of Dennison, Ginsberg and 
',7ebcr M he concluded that this value of Kr probably should bo increased 
thirty-fold. The water vapor radiative transfer problem of this thesis 
used the Blsasser Atmospheric Radiation Chart which, in turn, is based 
primarily on the new absorption coefficients mentioned above* Actually 
Kji was not computed in this study but the results of this study form a 
basis which would permit the calculation fairly simply. 

It was found also that carbon dioxide plays a surprisingly large part 
in the radiative heat tronsfor processes near the ground where rather large 
lapse rates are foxind. In the preparation of his radiation chart, Blsasser 
made the assumption that the carbon dioxide contribution to the not flux 
could be neglected in comparison to that duo to water vapor. It turned 
out in this study that the radiative contribution to the rate of change 
of temperature duo to carbon dioxide was in the neighborhood l/j to l/4 
that due to water vapor# Elsasser’s remark, concerning the effect of carbon 
dioxide, "that even thin layers radiate as blacl: bodies within these bands" 
was shotm to bo definitely invalid. In fact, using a formula for carbon 
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dioxide absorptivity due to Callendar [ 5 ] , the computed absorptivities 
ranged anywhere from 2 fo to 20^ for layers near the earth's surface, de- 
pending upon their thickness* 

Surprisingly largo values of turbulent temperature-change por half 
hour were obtained, in this study, at three inches above the ground* For 
example, there was a turbulent rate of cooling of 10°C per half-hour at 
1J00 local tine. It was found, however, that oven such large values wore 
negligible insofar as affecting the magnitude of the coefficient of eddy 
conductivity in the surface layer. In other words, whenever tho lapse 
rates wore significantly different from dry adiabatic, it was found that 
the turbulent tcnperature-change term, of the partial differential equation 
of oddy conduction 7/as negligible in comparison to the turbulent heat flux 
through either boundary of a slab of air* Thus the surface layer extended 
at least to 80 feet throughout the period of this investigation and during 
the night extended to about 200 feet* 

Other results deduood in this study are the following: The coefficient 

of oddy conductivity in the surface layer appears to follow a power law of 
the form K = Kq Z n , with n the exact value of m depending in sone way 
upon the stability, lloreover in this layer, K takes on its maximum value 
at about 1J00 LOT and its minimum value at about 1?00 LOT. 
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II. THE NATURE OF TEE DATA 



1. The Data. 

Data utilized in this research -was taken from Deport Ho. 29 com- 
piled by the Electrical Engineering Research Laboratory of the University 



and vapor pressure for selected levels from tho ground up to and including 
238 feet. The soundings are recorded at half-hourly intervals throughout 
the 24 hour period chosen to represent a particular synoptic situation. 
Ileasurements of wind speed and direction are recorded for selected levels 
up to 507 feet. The sounding site is on level ground under cultivation. 

The authors of M give the following statements of accuracy of the 



(a) Tho temperatures are accurate to within J. 0 C over tho 0-50°C range. 



sufficient to remove all turbulent fluctuations; however, they do not deviate 
more than 1.2°C from the >0-minute mean temperature. This turbulent fluctu- 
ation of the 6 -minute mean as compared with the JO-minute mean will be called 
"the gust error” , 

(c) The Aorovanes show a uniform response with an accuracy of 12/'. 

For a detailed description of the data and the field installation, tho 
reader is referred to the report of the data M. 

2. Selection of Data for this Study, 

The data for the period September 27-28, 1948 was selected for this 
study. A high pressure cell dominated the area and gave conditions of 
horizontal homogeniety. 



of Texas 




This report contains recorded measurements of temperature 



data. 



(b) The o-minute averaging period used to obtain tomperatures is not 
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Soundings at three-hourly intervals were used and computations of 
heat flux were made for the layers of air hounded by the ground, 6 and 
18 inches; 3, 6 , 12, 20, 35, 55, SO, 110, 145, 105, 235 and 288 feet, 
inasmuch as these are selected levels for which the temperature data 
were available. 
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. THEORETICAL A1ID EMPIRICAL CONCEPTS 

1, The Water Yapor Radiative Transfer Problcn. 

Elsasser [ 5 ] gives a method for determining the net flux due to 
the water vapor content of a layer of atnosphero through the use of his 
radiation chart. The corrected optical depth of water vapor, u, in a 
column of air of unit cross-section is computed by using the formula: 



l 

V. ~ ~ /coo 




(3.1) 



When only temperature and moisture data are available at given levels, 
prossures for these levels nay bo computed from the surface pressure by 
use of the integrated form of the hydrostatic equation: 



g AZ = Rj T (An. jj>, - f>a). (5-2) 

For the evaluation of net flux into or out of a slab of air, utilizing 
aroas on Elsasser’s Atmospheric Radiation Chart, Elsasser recommends the 
measurement of these aroas by means of a planineter* The difference of the 
net fluxes for a slab is then given in gram-calories per three hours by this 
method and i3 converted to rate of change of temperature of the slab by means 
of the equation: 
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In this equation, F 2 is the net upward flu:: at the top of the slab, while 
is the net upward flux at the bottom. of the slab. For further details 
concerning this method, reference should bo made to Hlsassor’s study. Heat 
Transfer by Infrared Radiation in the Atmosphere [5J • 

2. The Carbon Dioxide Radiative Transfer Problem. 

Callcndar £5] gives formulas for the absorptivity in terns of the 
quantity of carbon dioxide for various wave bands of a parallel beam of 
radiation. His formulas for the equivalent thickness, c, of carbon dioxide, 
and for the monochromatic absorptivity of carbon dioxide, in a column of air 
of thickness Z, pressure P & (in atmospheres) and temperature T (degrees 
absolute) are: 



c =8.7 Z f>i/T , M 

R-- I ~ + . ( 5 . 5 ) 

In equations (j> .4) and (?.5)> c is the length of absorbing path in crus of 
carbon dioxide at II.T.P. and is based on air having a normal carbon dioxide 
content (the partial pressure of carbon dioxide in the air now being .00p2 
atmospheres). The constants n and x are given for different limits of the 
intense absorption band centered at 15 microns. 

Elsasscr \_3~J states that, in the case of square-root absorption, a 
slab of thickness u is mathematically equivalent for diffuse radiation to 
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a linear column of length 1.7 8 u for beamed radiation. Tho condition of 
square-root absorption is satisfied for slabs of air of thickness 4 to 100 



should be increased slightly; but the error introduced for such layers Trill 
not bo great if the factor 1.7o is usod. Inasmuch as diffuse radiation 
occurs in the atmosphere, slab thicknesses must be increased by the factor 
1.78; rrhich amounts to an increase in the equivalent thickness, c, of 
equation (3.4) by the factor 1.78. Utilizing this corrected equivalent 
thickness in the absorption formula (3*5) allovrs this formula to be applied 
to slab thicknesses of the atmosphere. 

In this study, the rate of cliange of temperature of layers of air near 
the ground due to carbon dioxide as well as water vapor is to be determined. 
Furthermore, since the Atmospheric Radiation Chart of Elsasser is to be used 
for the rate of temperature-change resulting from water vapor, consideration 
must be given to his treatment of carbon dioxide* Elsasser [5] considers 
that the absorption by carbon dioxide occurs, for the most part, in the range 
752 cm”! to 584 cn"l, that is, between lp.p and 17.1 microns. He further 
considers that "very thin layers radiate as black bodies" and, therefore, 
neglects the radiative effects of carbon dioxide in comparison to those of 
•water vapor, for problems concerning net flux. The chart, however, does 
give a separate representation of the black body flux contained within the 
strong band interval 13«3 - 17 .1 microns. Callendar [ 5 ] gives the values 
n = .35 and x ■ .55 as suitable for determining the absorptivity by (5) 
within a band of limits 15*0 - 17.0 microns. It is assumed that the values 
of monochromatic absorptivity. A, determined from (3.4) and (3.5) are appro- 
priate for the carbon dioxide limits of Elsasser, 13 .3 - 17.1 microns. 



moters according to Panofsky 




For thinner layers, the factor 1.78 
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Knowing the ground temperature and the mean temperature of the 
layers of atmosphere under consideration, the flux through the top 
boundary of each layer can be computed* For example, for the first slab 
above the ground, the equation of radiative transfer, sometimes called 
Schwarzchild*s equation, gives 



E - (1-Al) V +A l E JU ' 



(5.6) 



mere E is the flux through the top boundary of the layer, Aq the mono- 
chromatic absorptivity of the layer, E^ 0 the black body flux in the band 
13*5 - 17*1 microns at the temperature, T, of the ground (or the flux into 
the lower boundary) and E^q is the black body flux at the moan temperature 
T of the layer. Strictly speaking Schwarzchild’s equation (3*6) is true 
for monochromatic intensity, or beamed flux, of radiation. Actually 
Elsasser £ 3 ] integrates the Schwarzchild equation in order to obtain tho 
flux emerging from an arbitrary slab. He finds that this integration in- 

CO ^ ^ ^ . 

volves the computation of integrals Ei 3 (x) = /p' e-’ dp where 

^xsec y/ f values of which are given in f5]. However, ho also shows that 
these computations imply the increase of tho absorbing path by the factor 1*78 
discussed above* It is therefore assumed that S chwarz child* s equation holds 
for fluxes E, E^ 0 , and E^q with the value of absorptivity A corresponding to 
tho absorbing path c replaced by 1*78 c, a result which Elsasser justifies* 

For slabs above the first, equation (3*6) also permits the computation 
of upward flux E through tho top boundary. However, for such computations 
C A= must bo taken to represent tho flux into the bottom of the slab, which 
was that out of the top of the slab immediately below, would still 
represent the black body flux, in the band 13*5 “ 17*1 microns, corresponding 
to the mean temperature of the slab under consideration. 
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In order to apply tho same procedure to the problom of computing 
fluxes for the downward radiation of carbon dioxide as for tho upward 
radiation, it is necessary to know the flux entering the top boundary of 
tho highest layer* The top boundary is 288 feet for this study* Through- 
out the period used in this study, whenever a ground inversion oxists, this 
height is very near tho top of the inversion. The slope of the sounding 
above the ground inversion during this period is quito similar to that of 
the standard atmosphere* This remark is, on tho whole, true also of the 
slope of the sounding above the 288 feet level during the hours in which 
no ground inversion is present. Therefore, tho atmosphere from 10 kilo- 
meters is divided into slabs of approximately ICO noters thickness and 
having the appropriate temperature according to the standard atmospheric 
lapse rate. In treating equation (5*6), slab 1 is to be considered the 
topmost slab, slab 2 that immediately below, etc* It will be shown below 
that the amount of flux from above 10 kilometers produces negligible 
effect at the 288 foot level. Let flux E Q entor the top of slab 1* Tho 
fluxes leaving the bottoms of the next several slabs are given by the 
following table. 

Slab Flux through the lower boundary 

1 ie o h, 

2 ml E fl + na 4 * b, (5*7) 

3 nml E + nma + nb c. 

o 

In (3.7 ), 1, m, n etc. correspond to tho factors. (l-A^), (l-iL,), 

(1-Aj), etc., respectively, of equation (3*6); a, b, c correspond to 

the values A-j_ E^, * ^3 *•*• e ’* :c * (3.6). 
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Since, using equations (J.4) and (J.5), the average value of 1, m, 

n, etc. is of the order 0.5 for a layer of thickness 100 meters, 

and E Q is of the order of 10 cal cm” 2 (5 hr.)'*' 1 ’, it turns out that the 
contribution of the initial flue: E Q is negligible in treating equation 
(J.6) for the downward directed flux. V/'e may therefore arbitrarily set 
S Q s 0 in determining the flux through the bottom of the slabs and still 
obtain the correct result at 288 feet. 

Finally, since we now may determine the amounts of upward and downward 
directed flux at every level, one can compute net-upward, net-downward and 
consequently net-outward directed flux, - F^, for every consecutive pair 
of levels. This determines the radiative heating, or cooling, of the slab 
due to carbon dioxide radiation, and the actual rate of temperature-change 
can be computed, as before, from equation (5*5). 

5* Heat Transfer by llolecular Conductivity. 

The upward flux of heat due to molecular conduction is given by the 
formula: 

jr 

F w ' ' * dZ ; (5-S) 

where oC is the coefficient of molecular conductivity of air, and has the 
value oC = 549(10)“^ cals cn“^ sec”^ T"^ at 0°C, but is a function of tem- 
perature (c.f. for example, p. 44 of M )• 

It is necessary' now only to compute upward fluxes F r 2 and F^ for the 
various slabs for which lapse rates arc known, and again use equation (J.?) 
to convert net flux out of a slab into a corresponding rate of change of 



(ID 



temperature. In actual practice it is normally found, since is very 
snail, that in order to register a measurable net flux, tho changes in 
lapse rates along tho vertical must be quite extreme, and this can be 
expected only within the lowest two or three slabs. 

4. Hate of Change of Temperature by Eddy Conductivity. 

The formula for the observed rate of temperature change is: 



2X , 


/IX) + 


/|X) W 


(AX) 




7>t 


l *r/w 


UWc v 


a* An VM/r, 


(3.?) 



where 




is due to inter vapor. 



is due to carbon dioxide. 



is due to molecular conduction 



> 




is duo to eddy conduction. 



The observed rate of change of temperature per Iialf-hour is given for 
each layer by the data used in this study. As a result of methods previously 
described, the only unknown of formula (3.0) is the rate of change of tempera- 
ture due to eddy conduction. 

5* Tho Coefficient of Eddy Conduction. 

The net-outward flux for a layer of atmosphere due to turbulent con- 
duction may be expressed as: 
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In (5.10) K2 and K-j_ are the coefficients of eddy conductivity at the top 
and bottom of any slab of atmosphere, respectively. Making the substitu- 



tions 



a -±- = J3 k. 

f RT ' RT 



and 



1©_ 



f 



(§f^) , 



both approximations being valid to Tilth in 1% for the layers between the 
surface and 288 feet, we obtain: 
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In (5.11), all items labeled with subscript 1 refer to measurements made 
at the center of slab 1, all items with subscript 2 refer to the center of 
slab 2; p^ is the mean of p-^ and p^. Thus (5.11) determines the rate of 
change of temperature approximately at the boundary of slab 1 and slab 2. 

In general, slab 1 will denote any lower slab and slab 2 the slab immed- 
iately above slab 1. Thus (J. 11) may be used as a recursion formula to 
determine all values of Iij. after has been determined. The only limi- 

lj > that is, when the lapse rate 
of tho slab is dry adiabatic, or, not significantly difforent from • 

In using (5.11) as a recursion formula, it is necessary to determine 
a value of K at some specific level to permit determination of all other Kj_. 
In order to do this without making too many assumptions, a value of K at 



tation on its use arises when - 
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some standard level can be determined by reference to the theory of 
momentum-transfer within the surface layer. The coefficient of eddy 
viscosity can be determined at a standard level, say 16 feet, and it can 
be assumed to give a representative value of the coefficient of eddy con- 
ductivity. The particular values of K vn.ll then depend to some extent on 
momentum-transfer theory, but the comparative values reflect only the 
measurements of eddy conductivity. Recent results of Lottau [>] are 
draTm upon for the purpose of finding a value of K at some standard level. 

According to the Prandtl theory of momentum-transfer*, ono obtains: 






( 3 . 12 ) 



^ U 



In (3.12) 1 is the mean mixing length and is the observed shear of 

the horizontal wind, both measured at a given level. Lettau r>] gives the 
following formulas which permit the calculation of the coefficient of eddy 
viscosity within the surface layer: 

, ^4 = ^( Z + Z °), <5*U) 



where 




(3.14) 



*Haurwitz, B. Dynamic Meteorology. Equation (74.2 ). LIcC-raw-Hill 1?41. 
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In (3,13) and (3*14), the subscript "a" indicates an isotropic or adiabatic 



3 e 

stirface layer, I n is the neon temperature of the slab for vrhich ^ 



been determined* The value of 
logarithmic law 






may bo determined by the well-knovm. 




In (3,15) and (3»l6), all variables are understood to refer to the surface 
layer. For example, 

r is the constant value of the stress, 

a 

u a is the horizontal wind speed, 

l a is the mixing length, 

k is the Von Ilarmnn constant, k - ,45, 



Z Q is the roughness parameter. 



The values of Z. and 



V 



t ; 



bu ; 



— can be obtained from a plot of 



vmcn 



adiabatic conditions are found. Observations of Lettau [ 8 ] indicate that 
such conditions exist at sunrise and sunset. An average u a profile, there- 
fore, caii bo drawn based on both sunrise and sunset data in order to detcr- 

(-a 



mine Z Q and the mean value of 



for the day. 
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IV. RATES OF CHANGE OF lENPERATUNE 



1 . Rate of Change of Temperature Due to T/'atcr Vapor • 

Radiation transfer due to vra.ter vapor was computed for each layer of 
the soundings chosen from the data "by the method described in Chapter III, 
Section 1, of this study. This method requires a knowledge of the specific 



level. These vapor pressures were converted to specific humidity by use 
of the Kiefer Isulti-Prossure Hygrometric Chart, Since equation (J.l) re- 
quires the pressure-thickness , £p, of all layers it -mas necessary to com- 



faco pressure must be known. The surface pressure was computed from the 
sea level pressure recorded for Austin, Texas on the surface map at the 
middle of the period under study, 27-28 September 1948. This was done by 
reversing the procedure of the reduction of pressures to sea level, knowing 
the height of the base of the tower (525 feet ISSL). 

A (u,T) relationship was computed for all levols of the soundings and 
plotted on the Elsas ser Chart. The difference of the net fluxes of each 
slab of the soundings was obtained by making tlio required area measurements 
with the aid of a plonimetcr. Figures 1 and 2 show schematically the areas 
that represent the differences of net flux, for one slab of atmosphere for 
tho following cases: 

(a) The slab lies within a ground inversion (Figure 1). 

(b) The temperature decreases with height (Figure 2). 

In Figures 1 and 2, tho reference levels and L ^ denote respectively the 
bottom and top of the slab, both taken to lie below tho 288 foot level. 



humidity at each level. Tho data 




vapor pressure at each 



pute the pressure at every level by equation (J.2). To do this, the sur- 
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Difference of Hot Fluoces for a Slab Belovr 288 Feet 

Figure 2. 
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Consider next Figure 1. ihe (u,T) curve extending above L-|_ must 
eventually cross that above L 2 because the former curve eventually corres- 
nonds to a greater integrated optical depth of water vapor. Suppose this 
"cross-over” point is laboled A. At 2200 LCT, the complete radiosonde was 
available to supplement the 2200 nicrometeorological sounding. It was found 
that the cross-over points for the (u,?) plots above the slab levels and 
L 2 occurred from about pOO feet for the lower, thinner slabs to about 800 
feet for the higher, thicker slabs. Threover, it was found that the area 
computation made by terminating the (u,T) curves above L-j. and L 2 at points 
3-j_ and B 2 corresponding to the 288 feot level was in error only by 2 - 1 % 
when compared with the areas given by the complete (u,T) curves* This small 
error is due to a tendency of positive and negative areas below the isotherm 
3]_B 2 to cancel each other. 

Consider next Figure 2, representing typical daytime conditions. Ho 
exact comparison is known here of tho error made in neglecting the area 
below and to the right of the isotherm. B^B 2 . Ho cross-over point occurs in 
this case so there is no tendency for cancellation of positive and negative 
areas. For the lower, thinner slabs little error could occur because the 
(u,l) curves above the slab boundaries are in general nearly coincident at 
and Bg. The greatest error should occur in the slab 185-255 feet and the 
error should be progressively smaller in the slabs immediately below. The 
neglected area, if appreciable, should cause a reduced amount of cooling or 
increased warning in the midday hours in the slab 185-255 feet. Actually 
Table 1 shown no evidence of reduced cooling of the slab 185-255 feet. For 
tho ratio of the 1500 LCT cooling of this slab to those of the threo slabs 
immediately below is almost identical with tho corresponding slab ratios of 
the 2200 LCT sounding, whose values arc known to be correct. 
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(Tines given in LCT . ) 
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(in feet) 
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Heat Transfer Due to A'ater Vapor 
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(Times given in LOT.) 



Layer 
(in feot) 


0200 

(°c) 


0500 

(°C) 


0800 

(°C) 


1100 

(°C) 


1500 

(°C) 


1600 

(°C) 


1900 
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2200 

(°C) 


0-.5 


-2.28 


2.28 


-6.85 


- 9.65 


- 5.47 


.00 


10.25 


.00 


. 5 - 1.5 
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2.55 


2.58 


1.51 
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- 1.24 


110-145 


- .45 


- .25 




1.80 


1.97 


1.04 


- .70 


-1.20 


145-185 


- .54 


- .70 


.55 


1.75 


1.65 


.81 


- .58 


- .85 
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- .65 


.41 


1.05 


1.14 


.70 
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Three-hourly Sate of Change of Temperature Due to hater Yapor 



TABLE 2 . 



( 21 ) 



\ - 



Table 1 gives the values obtained for heat transfer duo to water 
vapor in cal. (3 hr.)"^» Table 2 gives the three-hourly rate of change 
of temperature for each layer of each sounding, computed from Table 1 by 
use of equation (3»5)* 

2» Rate of Change of Temperature Duo to Carbon Dioxide. 

The differences in upwar d-dir acted flux due to the radiative transfer 
of carbon dioxide, in the band 13 *3 - 17 *1 microns , were computed for each 
layer by the use of equation (3*6). The differences in downward-directed 
flux were obtained from the method of Chapter III, Section 2. 

The computations of downward-directed flux took into account that the 
mid-period sounding vre is, on an average, ten degrees warmer throughout than 
that of the standard atmosphere, although the slope was assumed to be that 
of the standard atmosphere. This amounted to increasing the average emission 
from the layers by the factor 10/£, which factor is the percentage increase 
in flux of equation (3*6) due to the increase in the moan temperature, 

T n , of each layer compared to the standard atmospheric temperature. The 
adjusted value obtained for the flux through the top boundary of the layer 
288-235 feet was 18.75 cal cm“^ (3 hr.)“^. This value was used as a constant 
value, for the period of this study and computations of flux differences for 
the layers of each sounding under study were made using equation (3*6). The 
values of black body flux of emission for the carbon dioxide band 13*3 - 17*1 
microns, at the mean temperature of oach slab were taken from the top of the 
Atmospheric Radiation Chart. The values of monochromatic absorptivity wore 
computed using equation (3*4) and (3*5)* 
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(Times given in LOT.) 



Layer 0200 0500 0800 1100 1300 1600 1900 2200 



(in feet) 


A? 


AF 


A? 


4F 


AF 


A? 
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1 
* 

\n 
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Heat Transfer Dtae to Carbon Dioxide, Cal. (3 hr.)**' 1 
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.62 


.75 


l»66 


7.52 


7.50 


5*47 


2.42 


.55 


1. 5-5.0 


.11 


.17 


2.07 


7.81 


7.25 


5.45 


.52 


.57 


5.0-6 


-.14 


-.12 


1.57 


5 .84 


5.20 


2,41 


-.57 


-,06 


6-12 


-.28 


-.50 


.55 


4.10 


5.65 


1,60 


-.55 


-.25 


12-20 


-.17 


-.18 


.76 


5*15 


5.25 


1.51 


-1.07 


-.25 


20-55 


-.08 


- .06 


,6o 


2,08 


2,24 


.87 


-.87 


-.24 


55-55 


-.12 


-.06 


.54 


2.15 


1.66 


.67 


-.50 


-.55 


55-80 


-.10 


-.05 


.47 


1.08 


1.15 


.50 


-.78 


-.56 


80-110 


-.16 


-.05 


.25 


.80 


.55 


.55 


-.66 


-.71 


110-145 


-.50 


-.08 


.27 


.62 


.75 


.15 


-.62 


-.65 


145-185 


-.56 


-.14 


*55 


.56 


.45 


.00 


-.61 


-.61 


185-255 


-.55 


-.25 


.50 


.08 


.22 


-.15 


-.62 


-.57 


Three-hourly 


Rate of 


Temperature -change 


Due to 


C arbor. 


Dioxide 





TABLE 4. 
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Tabic 3 gives the results obtained for net-outward radiation. 

Table 4 rives the resultant rate of temperature-change for each slab 
of each sounding, computed using the results of Table 3 and equation 
3* Resultant Radiative Rate of Temperature-Change. 

The temperature-change values of Tables 2 and 4 were added algebraically 
yielding a resultant radiative tenperaturo-change for the slabs. Since 
these tenperature-changes are to bo compared with the observed tenperature- 
change at the various boundaries it was necessary to adjust them to these 
boundaries. This vras done by considering the radiative temperature-changes 
as existing at the midpoints of the slabs and then interpolating to the 
slab boundaries on log-log paper. These results are presented in Table 5* 
vrheroin all three-hourly rates have been reduced to half-hourly rates by 
dividing by 6. 
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(Times given in LOT) 



Boundary 
(in feet) 


0200 

(°C) 


0500 

(°C) 


0800 

(°c) 


1100 

(oc) 


1300 

(°c) 


’ 

■ 

St 


1300 

(°C) 


O'"' 
0 0 

OJ O 
CM — 


.5 


-.20 


.24 


-.40 


2.21 


2.77 


1.42 


2.02 


.03 


1.5 


-.77 


-.75 


1.15 


3,23 


7.68 


3.71 


•35 


-.58 


3*0 


-.71 


-.65 


1,43 


6.13 


6.34 


2.33 


-.32 


-.51 


6.0 


-.47 


-.47 


1.40 


4.25 


3.72 


1.38 


-.42 


-.53 


12 


-.50 


-.55 


.75 


2.88 


2,60 


1.35 


-.57 


-.27 


20 


-.15 


-.24 


,47 


1,87 


2.06 


.35 


-.57 


-.20 


55 


-.12 


-.12 


.34 


1,28 


1.25 


.52 


-.45 


-.15 


55 


-.06 


-.05 


.51 


.35 


.86 


.40 


-.41 


-.13 


8o 


-.04 


-.06 


.23 


.60 


.53 


.30 


-.51 


-.28 


110 


-.08 


-.05 


.13 


.47 


.50 


.24 


-•26 


-.32 


145 


-.15 


-.10 


.14 


.38 


.41 


.17 


-.21 


-.28 


185 


-.15 


-.14 


.17 


.28 


.30 


.12 


-.20 


-.23 


255 


-.15 


-.11 


.10 ' 


• 16 


.20 


.06 


-.18 


-.13 



Half-hourly Rato of Temperature-change: Water Vapor + Carbon Dioxide 



TABIJ3 5. 
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4, Hate of Change of Temperature Due to Holecular Conduction* 

The flux due to molecular conduction vras computed by the use of 
equation (2*8). The temperature, riven in the observed data, for each 
level was used, hence the flux was considered as passing through the 
center of each layer. Consequently, the temperature-change computed 
from the value of the flux differences for two consecutive layers must 
be considered as being at the boundary between the two layers. This 
resulted in a rate of temperature-change at the boundary of each slab 
without the necessity of interpolation. The results obtained for rate 
of temperature -change due to molecular conduction are given in Table 6, 

(Times given in I,CT) 



Eoundary 
(in feet) 


02 C 0 

(°c) 


0500 

(°C) 


0800 

(°0) 


1100 

(°C) 


1300 

(°c) 


1600 

(°r) 


1900 

(°c) 


2200 

(°C) 


• 5 


-.?o 


- 1.90 


-i.6o 


7.20 


8.0 


4.10 


-.57 


-2.60 


LT\ 

• 

I — 1 


-.12 


-.16 


-.09 


.97 


.45 


.17 


-.27 


0.0 


3.0 








»o6 


.07 


.05 


-.06 





Half-hourly Rate of Temperature-change Due to 
Holecular Conduction 

TABLE 6. 

5 . Rate of Change of Temperature Due to Eddy Conduction, 

Those results follow imnodiately by the use of the observed 
and the values of Table 5 and 6 in equation (>?). The values for the 
rate of temperature-change due to eddy conduction are given in Table 7. 



( 27 ) 



(Tines given in LCI.) 



Boundary 
(in feet) 


0200 

(°c) 


0500 

(°C) 


0800 

(°c). 


1100 

(°c) 


1300 

(°c) 


1600 

(°c) 


1900 

(°c) 


2200 

(°c) 


.5 


. 8 o 


1.66 


l. 6 o 


-3,01 


-10.37 


-6,32 


-1,05 


-1,67 


1.5 


0.55 


0.75 


1.22 


- 7 . 6 o 


-7.63 


-4.48 


-.02 


0.18 


5 


0.21 


0.45 


0.47 


-4.95 


-5.91 


-5 .48 


-.02 


0.01 


6 


-0.2J 


0.17 


0.50 


-3.15 


-3.62 


-2.28 


-0,28 


-0.07 


12 


-0.30 


0.05 


0.57 


-2,10 


-2.20 


-1.23 


-0.33 


-0.23 


20 


-c.31 


0.04 


1.43 


-1.57 


-2,36 


-1.15 


-0.03 


-0.30 


35 


-0.48 


-0.28 


1.46 


-0.78 


-0.85 


-0.52 


0.15 


-0.55 


55 


-0.34 


-0.11 


1.4? 


-0.45 


-0.56 


-0.30 


0.31 


-0,21 


8o 


-0.06 


-0.04 


1.57 


0.00 


-0.5? 


-0,10 


0.21 


-0,22 


no 


0.28 


0.05 


0.47 


-0.07 


-0.10 


-0,14 


0.06 


-0.28 


145 


0.03 


0.30 


1.46 


0.02 


-0.01 


0.03 


0.01 


0.28 


185 


-0.05 


0.14 


1*33 


0.02 


0.20 


0.08 


-0.10 


0,43 


235 


-0.35 


0.12 


1.30 


.04 


0.10 


0,04 


0.08 


. 3 ? 



Half-hourly Tenperatur e-change Due to 
Dddy Conduction 



TABLE 7, 
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V. RESULTS ALTD CORCLUSIOIIS 



1, I'ethod of Computation of the Coefficient of Eddy Conductivity, 

As indicated in Chapter III, Section it was decided to determine 
for all soundings evaluated a first value of K from the theory of turbu- 
lent transfer of momentum. For this purpose it was necessary to make use 

du 

of observed winds in. order to evaluate dZ for use in equation (5,12), 
Two winds were given consistently in tho data which were always within the 
surface layer, one at 12 feet and one at 41 feet. For each of the eight 
soundings studied, these twx> values of the wind were plotted against log Z 
on seni-log caper, and a straight line joining tho two points was drawn, 

m 

The value of oZ in the vicinity of 16 feet vra.s then read off as being 

Au 



the value AZ corresponding to a height increment of 100 centimeters 
at 16 feet. The reasons for choosing the 16 foot level as tho standard 
reference level are two-fold: 

(a) The level had to be within the lowest levels for which wind 
data are available, 

(b) The level had to be within the lowest slab consistent with (a) 
for which lapse rate data is available. This is necessary since, in 

£©. 4 /ax . 

equation the nixing length, 1, involves * (^3% ay * 

In satisfying requirement (b), it was necessary to use the slab 12-20 

6 e 



feet, which was considered to define dZ at l6 feet. 

In several cases however, it ms necessary to shift to the level 2>«5 
feet. This necessity arose when an obviously inconsistent lapse rate ’was 
indicated in the layer 12-20 feet. This shifting to the level 2J.5 feet 
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amounted to dealing with the over-all lapse rate in the layer 12-35 feet, 
which then ;ave A7 at 25*5 feet. 

In order to evaluate the mining length, 1 (the remaining ±©rm. of 
equation (5,12)) at the standard level, a knowledge of the corresponding 
value l a , the mixing length in an isotropic atmosphere, is necessary. 
According to Lettau r»] isotropic or adiabatic turbulence is most nearly 
realized around sunset or sunrise. Thus, in order to derive the adiabatic 
velocity-profile for the 24-hour period under investigation, an average of 
three half-hourly wind-speed reports at both sunrise and sunset was obtained* 
Then, plotting mean wind-speed against log Z on semi-log paper for the two 
values of Z a straight line was obtained. Figure J>, which according to 
Lettau should represent the wind profile for Manor, Texas under adiabatic 
conditions for 27-28 September 1948, The straight line obtained apparently 
does ropresent a close approximation to the true adiabatic atmosphere, since 
Figure 3 indicates a roughness parameter Z - 2 ems, which checks rather well 
with values quoted by Stewart [ ll for fallow fields* Knowing the value of 
the roughness aerometer uermits the computation of 1 from ecuation ( 5 * 13 )* 

iv* 

The value of at lb feot was then read directly off the velocity 

profile of Figure 3» Since 16 feet ± 490 ems this can be done quite simply 
by reading off the wind-speed differences between bOO and 400 ems and divid- 
ing by 2C0 cns. This calculation leads to *15 sec*"^, a value which may be 
checked mathematically using the theoretical equation of the straight line 
of Figure 3» This is done below. 
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Figure >. 
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The equation of the wind profile line of figure 3 be written 

as in (3*15 )• Ur. in£ the nean value of U a at the 12 foot level (366 ems), 

-1 



1 



45, Z Q ■ 2 ems, U a ■ 3*8 tips, vre obtain - ?2 



♦ 8 



cm sec 



Then from equation ( 16 ), 



(■ 



but 

bZ 






,45 (52. ) ( 4 ^—) = .15 seo ” 1 . 



but 



Similarly at 23«5 feet the value of ^2 710,8 bo be: 

d-u& 



az 



= 743 *^ 2 . 8 )^^- = ,10 sec" . 



At these two levels l a = 220 and 323 ems, respectively. Values of K at 

de 

the reference level nay now be determined once the values of dZ » ©c 

6 u 

are kno-Rjfc. These results ere tabulated below in Table 8 for the 
various times of days 





0200* 


0 

0 

0 

* 


1100* 


1300 


1600 


1900 


2200 


M 

dZ 


m y &* , .oood 7 


,00038 


-.00073 


-.0026 


-.00061 


.00092 


,00092 


©c 


a< /cm *0050 


.0050 


.0031 


.0066 


0 

n\ 

O 

O 

• 


.0066 


.0063 


AiL 

dZ 


sec'* *183 


.34 


.213 


• 

-£* 

cx> 


.222 


.214 


CO 

N'n 

■'d" 

* 


K 


c^ec 11 ^ 000 


28,000 


40,900 


63,800 


38,000 


8,000 


16,000 




C omputat ion 


Values 


for Determination of K 


at Reference Level 





TABLE 8. 
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In Table 8, the * alongside the tine indicates that the reference 
level iras talcen at 2J.5 feet due to the non-representative lapse rates 
between 12 and 20 feet* The 0800 values of K were not computed since 
many of the slabs hal either non-representative cr dry adiabatic* lapse 
rates so that the heat transport equations (>,11) yield no solutions, 

For all soundings, other than the 0800 sounding, it was possible to 
obtain solutions of equations (3*11) usin'-; the values of K of Table 8 as 
a starting point. The values of h for 0200, 05CO, 1100, lJOO, l60C, 19C0 
and 2200 LOT at the indicated elevations are listed in Table ?, and shown 
graphically in Plate I, In Table 9 the occurrence of blarf: spaces indi- 
cates the existence of a non-representative lapse rate in the slab having 
the indicated height as its midpoint. It was therefore necessary to con- 
solidate adjacent slabs, obtaining a value of K corresponding to the new 
midpoint, as discussed earlier in Section 1 of this chapter, 

♦Actually, isothermal lapse rates vrere not significantly different from 
the dry adiabatic for many of the layers under study, when the gust error 
of ± .2°C is considered. 
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(Times given in LOT.) 



Level 

(feet) 


0200 

L ( 10)2 


0500 

X (10 )2 


1100 

K(10)3 


1300, 

r(io)2 


1600 

k(10)3 


1900 

k(io)5 


2200 

Yv(ioy 


.25 


.075 


.14 


.045 


.26 


.04 


.036 


.099 


1.0 


.51 


.68 


.50 


3.7 


.66 


• l£ 




2.25 


1.9 


5.0 


2.4 


16.4 


1.7 


.41 


2.3 


4.5 


1.5 


2.4 


c - 0 


180. 




•73 


2.3 


7.5 










15. 






5.0 


4.9 


5.9 


20 . 


48. 




1.3 


8.7 


16.0 








64 . 


38 . 


C.o 


16. 


25.5 


15. 


28. 


41. 










27.5 








500. 


63 . 


5.1 


?.2 


45.0 


11. 


2 ?. 


480, 


1700. 


35. 


8.0 


11. 


67.5 


17. 


28. 




300. 


130. 


32. 


9.0 


P2.5 






120. 










95.0 


8 . 


21. 








73. 


20. 


127.5 


10. 


18. 


550. 






25. 


530. 


165.0 


17. 


20. 








28. 


840. 


210.0 


17. 


62. 








210. 


920. 


261.5 


10. 


87. 








190. 


1500. 



The Coefficient of Eddy Conduction* K, (cm 2 sec“") 



TABLE 9. 
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In order to show moro clearly the variations of K with both height 
and time, graphs of K at fixed elevations were drawn showing the tine- 
variation of K at these elevations. The fixed elevations used are 3 
inches, 6, JCj 60 and 125 feet and the results are shown in Figures 4, 

5 , 6 , 7 , and 8 . 
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Diurnal Variation of "K" 



at .25 Feet 



Figure 
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Diurnal Variation of "K" at 



30 Feet 



Figure 6. 




0 <xo 
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Diurnal Variation of "K" at 60 Feet 



Figure 7 . 
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Diurnal Variation of "K" at 125 Feet 



Figuro 0. 
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2. Discussion of Results, 

Lettau [?] has defined the surface layer essentially as the layer 
in which the turbulent time-rate of change of a property is negligible 
in comparison with its turbulent flux. In the solution of equations ( 5 . 11) 
for the various levels, it is evident that the value of the lapse rates be- 
come exceedingly critical as the value of the lapse rate approaches the dry 
adiabatic . For, if the maximum gust error ^",2 °C occurred in opposite 
senses at adjacent levels, it would cause an error in observation of the 
lapse rate by sevoral tines ^ , depending upon the layer thickness. Thus 
it follows that the most representative values of K may be expected where 
the lapse rate If is largest compared to and especially where it is of 
the order of 50-100 (or more) tines . In order to emphasize this de- 
pendence upon lapse rate, let us consider the values of the lapse rate in 
the various layers for the 0200 and 1^00 LOT soundings. These are shown 
in Table 10. In this table, the abbreviation "n.r." indicates a non- 
representative lapse rate; therefore, necessitating a consolidation of 
the two adjacent layers, 

Ilotice, from Table 10, that at 130 0 LCT the lapse rates were strongly 
super-adiabatic near the ground but not significantly different from the 
dry adiabatic above 80 feot. At 0200 LCT the lapse rates were more nearly 
uniform than at IJOO LCT, A comparison shows that the lJOO LCT lapse rates 
were numerically greater than those for 0200 LCT between the surface and 6 
feet, but, above G feet, the early morning lapse rates in general exceed 
those of nid-day. " ? or purposes of comparison, the value of the dry adiabatic, 
lo.pse rate is = .98 x 10“ " °C per cm. 
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Lapse Rate (°c/cn) x 10“4 



Layer 
(feet ) 


0200 LOT 


1500 LCT 


0-.5 


-1.1x10? 


7.1x10? 


,5-1.5 


-1.6;:10 2 


4 • 5x1 0 2 


1.5-5 


-44 


110 


5-8 


-55 


11 


6-12 


-16 


58 


12-20 


n.r . 


25 


12-55 


-5.7 


15 


20-55 


n.r. 


4.4 


55-55 


-6.6 


4.5 


55-80 


-5.5 


1.5 


80-110 


1 

VO 

A 

CO 


8.0 


110-145 


-7.5 


isotherm* 


145-185 


-4.1 


II 


185-255 


-5.5 


It 



Values of Lapse Rates at 0200 and 1JC0 LOT 



TABLE 10. 
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In general, it was found that when Y was significantly different 
from , the turbulent rate of change term of equations (Jill) was 
completely negligible-, Due to the inability to smooth out the turbulent 
fluctuations of temperature, the definition of the surface layer in this 
study is the layer in which the lapse rate is significantly different from 
dry adiabatic, both the gust and instrumental errors being considered* Thus 
the surface layer at 1J00 LCT extends to about 80 feet, but at 0200 LCT ex- 
tends to at least 235 feet. This is evident from the values of Table 10. 

It must be realized that a negative value of this lapse rate is more signifi- 
cantly different from than a corresponding positive value of the lapse 

rate. Then, for all practical purposes, equation (3*11) within the surface 
layer reduces to: 



where subscripts 1 and 2 refer to a lower and upper level, respectively. 

Let us consider tho actual values of K shown in Plate I in which the 
value of K is plotted on a’ log-log scale against the height of the midpoint 



sounding, the first three points in order of ascending height are very close 
to a straight line, there being some deviation from tho straight line re- 
lationship as the fourth, fifth, etc., poinbs are considered* This deviation, 
commencing with the fourth point, is due mainly to two factors: 

(1) Decreasing, therefore less representative, values of the lapse rate* 

(2) The change in temperature elements above the three foot level; there 
being non-aerated resistors bolow and aerated resistors abovo. 




o 




of the slab to which the value of K corresponds. It is noted that, for any 
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Therefore, it was considered remissible to fit the best lire to each 
sounding in Plate I up to the maximum height for which II and Z both in- 
crease together. The line of best fit for each sounding was drawn in 
accor lance with the following principles: 

(a) Through the first point, 

(b) Through as many additional points as possible. 

(c) Separating approximately equal numbers of points. 

The slopes and heights to which these lines nay be considered representa- 
tive are shovm in Table 11 for each sounding. 



Sounding 
(Tine in LCT ) 



Slone of line, m 



T 'axinun height 
of 15:near distri- 
but ion (feet) 



Values of Kq 
in equation (5o) 
(erf " sec~^- ) 



0200 


1.14 


23.5 


365 


0500 


1,17 


23,5 


506 


1100 


1.79 


45.0 


536 


1500 


1.69 


45.0 


26^0 


1600 


1.70 


23.5 


410 


1900 


1.05 


45.0 


155 


2200 


1,25 


16,0 


544 


Slopes 


of Log K Against 


Log Z for the Lines 


of Plate I 



TABLE 11 
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There is a -tendency for the lines for the 1^00 and 2200 LCT soundings 
to continue to be representative up to 235 feet* Hoy/ ever , for the 0200 and 
0500 LCT soundings beyond 25 .5 feet, the points fall well below the line. 

This suggests the possibility of fitting a curve other than a straight line 
(perhaps a parabola) in tho early morning hours; however, this idea was not 
pursued further. 

If we write the equation of any straight lino of Plate I in the form 

log K = m(log Z) + b, (5.2) 

where n is the measured slope (see Table 11) of the line, we obtain upon 
raising to povrers of 10, the equation 

K - I0 b Z n = Kq Z m . (5.5) 

Since each line passes through the point corresponding to Z - 1/4 foot, the 
value of Kp for each sounding may be assigned by setting Z - 1/4 and K to 
its appropriate value. The values of IZq for each sounding are shown in the 
last column of Table 11. Equation (5.5), with the appropriate values of K] 
and m, represent the distribution of K with respect to height for the various 
times of day and up to the heights indicated in Table 11. This is the chief 
empirical result of this study. 

Finally, we shall state the conclusions resulting from the foregoing 
analysis and from Figures 4 through 8. These conclusions are: 

(1) The average mid-day value of n is 1.73. 

(2) The average value of ra near midnight is 1.18* 

(3) The value of m seems to approach unity at sunset, 

(4) The exponent m appears to be a function of stability. 
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(5) The maximum value of K, level for level up to 60 feet, appears 
to occur at about 1J00 LCT , at about the time of the maximum surface tem- 
perature, The power lav; satisfied by K at this time of day is 

K = 2660 Z 1 *^, 

where Z is the height in feet, 

( 6 ) The minimum value for II, level for level up to 60 feet, appears 

to occur about lJOO LCT, with a very slew rate of increase during the night, 

(7) From Figure 8 at 125 foot there appears to be a second pronounced 
diurnal maximum of K at 2200 LCT, This appears to be due to the diurnal 
variation in the height of the surface layer. According to Lcttau M. 
the surface layer is characterized by increasing values of K with height, 
with decreasing values of K occurring above this layer. If, during the day- 
light hours the 125 foot level is above the surface layer, but passes into 
it ‘at about 2200 LCT, the large increase in the value of V at 125 feet ap- 



pears to be reasonable 
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